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ABSTRACT 

The TeV blazar PKS 2155—304 was monitored with the X-ray sateUite 
ASCA in 1994 May, as part of a multiwavelength campaign from the radio 
to X-ray bands. At the beginning of the two-day continuous observation, we 
detected a large flare, where the 2 — 10 keV flux changed by a factor of 2 on 
a time scale of 3 x 10^ sec. During the flare, the increase in the hard X-ray 
flux clearly preceded that observed in the soft X-rays, the spectral evolution 
tracking a 'clockwise loop' in the flux versus photon index plane. Ascribing the 
energy-dependent variability to differential synchrotron cooling of relativistic 
electrons, we estimate the magnetic field B in the emission region. We tested 
two different methods of comparing the time series in various X-ray bands: (i) 
fitting the hght curves to a Gaussian and searching for the time shift of the 
peak of the flare, and (ii) calculating the discrete correlation function. Both 
methods yielded a consistent solution of 5 ~ 0.1 G. We also found that the flare 
amplitude becomes larger as the photon energy increases, while the duration 
of the flare stays roughly constant through the ASCA energy band (0.7 — 7.5 
keV). In the framework of the time dependent synchrotron self-Compton model 
in a homogeneous region, we consider a flare where the maximum Lorentz factor 
(7max) of the injected electrons increases uniformly throughout the emission 
volume. The temporal evolution of spectra as well as the hght curves were 
reproduced with the physical parameters self-consistently determined from 
seven observables. We obtain S ~ 0.1 — 0.2 G and a region size R ~ 10~^ pc, 
for relativistic beaming with a Doppler factor of 5 ~ 20 — 30. We discuss the 
significance of light travel time effects. 

Subject headings: keyword: BL Lacertae objects: individual (PKS 2155—304) — 
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X-rays: galaxies — radiation mechanism: non-thermal — gamma-rays: theory 
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1. Introduction 

Recent observations with the EGRET instrument (30 MeV - 30 GeV; Thompson et 
al. 1993) onboard the Compton Gamma-Ray Observatory (CGRO) reveal that more than 
60 blazars are bright 7-ray emitters (e.g., Mukherjee et al. 1997). Multi- frequency studies 
of blazars show that the overall spectra (plotted as uF,y) have at least two pronounced 
continuum components: one peaking between IR and X-rays, and another in the 7-ray 
regime (e.g., von Montigny et al. 1995). The strong polarization observed in the radio 
and optical bands implies that the lower energy component is most likely produced by 
synchrotron radiation of relativistic electrons in magnetic fields, while inverse- Compton 
scattering by the same electrons is believed to be the dominant process responsible for the 
high energy 7-ray emission (e.g., Ulrich, Maraschi, & Urry 1997; Ghisellini et al. 1998; 
Kubo et al. 1998). However, the source of the 'seed' photons for the Compton process 
is a matter of debate. In some models, it is synchrotron radiation internal to the jet, as 
in the Synchrotron-Self-Compton (SSC) model (Jones ct al. 1974; Marscher 1980; Konigl 
1981; Ghisellini & Maraschi 1989; Marscher & Travis 1996), while in others, it is radiation 
external to the jet, as in the External Radiation Compton (ERG) model (Dermer & 
Schlickeiser 1993; Sikora, Begelman, & Rees 1994). 

Observations with ground-based Cherenkov telescopes indicate 7-ray emission 
extending up to TeV energies for 4 nearby blazars: Mrk 421 (z = 0.031; Punch et al. 1992), 
Mrk 501 {z = 0.034; Quinn et al. 1996), lES 2344+514 {z = 0.044; Catanese et al. 1998) 
and PKS 2155—304 {z — 0.117; Chadwick et al. 1999). These extreme blazars are similar to 
each other in that their lower energy components peak in the X-ray band, while the higher 
energy peaks are in the TeV bands. Various authors have interpreted the multi-frequency 
spectra of TeV blazars by a simple homogeneous SSC model (e.g., Mastichiadis & Kirk 
1997; Pian et al. 1998; Kataoka et al. 1999a). 
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Rapid and large amplitude variability is a marked feature of blazars. Flux variations on 
a time scale of hours to days indicate that blazar emission comes from a very compact region, 
and is likely Doppler-boosted towards us as a result of relativistic motion of the emitting 
blob of plasma (e.g., Mattox et al. 1993; 1997). The correlation of variability between 
different bands can potentially discriminate between various models. The multi-frequency 
campaign of Mrk 421, conducted in 1994 (Macomb et al. 1995; 1996; Takahashi et al. 1996) 
revealed correlated activity between the keV X-ray and TeV 7-ray emission. Similarly, 
Mrk 501 showed a remarkable flare in 1997 April, where both X-ray and TeV 7-ray flux 
increased by more than one order of magnitude from the quiescent level. During the flare, 
the synchrotron peak was observed around 100 keV, which is a shift by a factor of more than 
100 from its location in the quiescent level (Catanese et al. 1997; Pian et al. 1998). For 
PKS 2155—304, TeV emission was detected only recently (Chadwick et al. 1999). Previous 
multi-frequency campaigns, conducted in 1991 November (Edelson et al. 1995) and 1994 
May (Urry et al. 1997), detected flares from UV to X-ray energies, with very different 
multi-frequency characteristics between the two epochs; there were no simultaneous TeV 
measurements. 

From the results of those multi-frequency campaigns, many authors have attempted to 
derive the physical parameters associated with the blazar emission and/or the origin of the 
flares. For example, Takahashi et al. (1996) associated the soft lag observed in Mrk 421 
during the 1994 observation with the energy dependence of the synchrotron cooling time, 
deriving a magnetic held strength S ~ 0.2 G for S — 5. Importantly, this was the flrst case 
where the magnetic field was calculated only from the observed X-ray spectral variability. 
However, it is not certain that the observed time lag can be straightforwardly associated 
with the the synchrotron cooling time. 

Mastichiadis & Kirk (1997) considered time-dependent SSC models where the 
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variability is on time scales longer than R/c, R being the size of the emission region. 
However, it is to be noted that the cooling time of the highest energy electrons, tcooi(7max), 
may well be shorter than R/c] in this case, rapid local variability could be relaxed 
(smoothed) by hght travel time effects over the source on a hght-crossing time scale. These 
conditions may be implied by the observed quasi-symmetric shape of the light curves seen 
both in X-rays and at lower energies (e.g., Urry ct al. 1997; Giommi ct al. 1998). Chiaberge 
& Ghisellini (1999) took this relaxation effect into account and applied their model to the 
Mrk 421 flare in 1994. 

In this paper, we independently develop a time-dependent SSC model, to follow the 
spectral evolution and variability patterns of blazars. In contrast to most previous works, 
our calculations are quantitatively compared with the observational data. We apply this 
model to the X-ray flare of PKS 2155—304 during the 1994 May campaign. Seven physical 
parameters for the one-zone homogeneous SSC model are self-consistently determined from 
seven observables and cross-checked with the magnetic field, derived directly from the 'soft 
lag' observed in the X-ray band. 

We present the ASCA observation and analysis of PKS 2155-304 in 1994 May in § 

2. Analyzing the time-series by various methods, the magnetic field is deduced under the 
assumption that the energy-dependent variability is caused by synchrotron cooling. In § 

3, we develop and summarize our time dependent model, and apply it to the fiare in 1994 
May. In § 4, we discuss the origin of the X-ray fiare and its relation to that observed at 
lower frequencies. Finally in § 5, we present our conclusions. 



2. Observations and Data Analysis 
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2.1. Spectral Fitting 

We observed PKS 2155-304 with ASCA during 1994 May 19.2 - May 21.3 UT, 
yielding a net exposure time of ~ 80 ksec (Makino et al. 1996). The observation was 
performed in a normal 1 CCD mode for the Solid-state Imaging Spectrometer (SIS: Burke et 
al. 1991; Yamashita et al. 1997), and a normal PH mode for the Gas Imaging Spectrometer 
(GIS: Ohashi et al. 1996). We applied standard screening procedures to the data and 
extracted source counts from a circular region centered on the target with a radius of 3 and 
6 arcmin for SIS and GIS, respectively. Since the count rate of the background (~ 0.01 
cts/s) and its fluctuation are neghgible compared with the count rate from PKS 2155—304 
in the 0.7 — 7.5 keV band, we performed no background subtraction for the data, to avoid 
any instrumental artifacts. 

The light curve obtained from the SIS detectors is shown in Figure 1. The data reveal a 
large flare at the beginning, followed by lower amplitude fluctuations. This flare includes the 
rising and decaying phase; both have a time scale of ~ 3 x 10^ sec. The source variabihty is 
somewhat different in different energy bands, and this is illustrated by separately plotting 
the light curves in the lower (0.5 — 1.5 keV) and higher (1.5 — 7.5 keV) energy X-ray 
bands. Notably, the amplitude of flux change is larger at higher photon energies - a factor 
of 2 at 1.5 — 7.5 keV, while it is a factor of 1.5 in the 0.5 — 1.5 keV band - but the flare 
duration seems to be nearly the same. Also note that the peak of the light curve in the 
hard X-ray band leads that in the soft X-ray band by ~ 4 ksec (Makino et al. 1998). 

The bottom panel of Figure 1 shows the time history of the hardness ratio, deflned as 
the ratio of the SIS count rates at 1.5 — 7.5 keV to those at 0.7 — 1.5 keV. The hardness 
ratio decreases as the flux decreases, indicating that the X-ray spectrum tends to steepen 
as the source becomes fainter. This is a general feature of HBL (high-frequency peaked BL 
Lac objects) spectral variability (e.g., Giommi et al. 1990). 
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To investigate the rapid variability of PKS 2155—304, we divided the total exposure 
into 5 ks intervals. For each segment, we analyzed the combined SIS/GIS spectra. A 
single power law function and the photoelectric absorbing column A^h fixed at the Galactic 
value (1.77 x 10^° cm~^; Stark et al. 1992), often used to represent X-ray spectra of 
BL Lac objects, do not fit any of the spectra well. A power law fitting the data at ;^ 2 
keV was too steep for the data at 2 keV. We found that a single power law plus a free 
neutral absorbing column is an acceptable model, with reduced ranging from 0.93/dof 
to 1.17/dof for ~ 1000 dof. We are aware that this is an unphysical model, since the source 
is a strong EUV emitter and the absorption cannot be as large as the ASCA data imply. 
Nonetheless this model is a convenient representation of the curved ASCA spectrum. 

We plot the best-fit spectral parameters for the 36 data segments in Figure 2. Figure 2 
(a) shows the time evolution of the differential photon index in 0.7 — 7.5 keV band, while 
Figure 2 (b) shows the change in the 2 — 10 keV flux. The spectral shape varied such that 
the differential photon index increased from F = 2.43 ± 0.02 to 3.07 ± 0.03. The flux also 
changed dramatically, with (14.5 ± 0.16) xl0~^^ erg cm~^ s~^ at the peak of the flare and 
(3.65 ± 0.06) xlO~^^ erg cm~^ s~^ near the end. Although we have allowed the 'equivalent' 
photoelectric absorbing column A^h to be free in the model flts, it stayed nearly constant at 
~ 1.0 X 10^^ cm~^ throughout the observation. 

This spectral change is best illustrated as a correlation between the flux and the photon 
index in Figure 3, showing that the spectrum is generally harder when brighter and steeper 
when fainter. Moreover the variations show a hysteresis, a 'clockwise loop,' clearly seen 
for the flaring period, followed by gradual steepening in the subsequent decay phase. This 
clockwise motion in the flux versus photon index plane has been seen in this particular BL 
Lac object during earlier X-ray observations (Scmbay ct al. 1993), and similar behavior 
was clearly seen during observations of Mrk 421 in 1994 (Takahashi et al. 1996) and H 
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0323+022 in 1987 (Kohmura et al. 1994). 



2.2. Correlation of the Variability in Various X— ray Bands 

Below, we quantify the energy dependence of the PKS 2155—304 X-ray light curves, 
concentrating on the intense flare from 1994 May 19.2 - 20.0 UT (indicated by arrows in 
Figure 1). 



2.2.1. Gaussian Fit of the Light Curves 

As Figure 1 shows, the X-ray hght curve during the flare of PKS 2155—304 is very 
symmetric. This implies that a simple technique can be applied to measure the energy 
dependent behavior of the light curves. We fitted the light curves in various X-ray energy 
bands with a Gaussian plus constant offset function. In this model, the photon count rate 
at an arbitrary time t is expressed as: f{t) — Cq + Ci x exp {—{t — ip)^/2(T^), where Co is 
a constant offset, Ci is the amplitude of the flare component, tp is the peaking time, and a 
is the duration of the flare, respectively. 

The light curves, binned at 1024 sec, were divided into 10 logarithmic-equal energy 
bands from 0.5 keV to 7.5 keV (for the SIS data) and 0.7 keV to 7.5 keV (for the GIS 
data; 9 bands). A Gaussian provides a sufficient fit of the resultant energy- binned light 
curves. In Figure 4, we show examples of the Gaussian fit, which turned out to be a 
good representation of the data; the probability of the fit was P(x^) > 0.1 for 12 of 19 
light curves. Note the difference in the peaking time and the difference of the ratio of the 
normalization (Ci) to the constant offset (Co). 

Figure 5 (a) shows the dependence on energy of the fiare duration derived from fitting 
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the energy- binned light curves. This result is consistent with a constant fit of cr = 1.5 x 10^ 
sec (sohd fine) . Figure 5 (b) shows the energy dependence of the flare amplitude, deflned 
as the ratio of the normalization of the flare to the constant offset: Ap = Ci/Cq. The 
flare amplitude becomes larger as the photon energy increases. This is mostly due to the 
decrease in the constant offset (Co) at higher energies. In other words, the spectrum during 
the flare is harder than that in the quiescent state (see Figure 2 (a)). 

In the subsequent analysis, in order to reduce errors and present the time lag more 
clearly, we binned the light curve into 5 energy bands. Figure 5 (c) is the lag of the peak 
time, calculated from the difference of of the Gaussian, as compared to that measured in 
the 3.0 — 7.5 keV band. As shown in the flgure, the hard X-ray (3.0 — 7.5 keV) variability 
leads the soft X-rays (0.5 — 1.0 keV) by ~ 4 ksec. We assume that the delay of the response 
of the soft X-ray flux reflects the synchrotron lifetime (cooling) of the relativistic electrons: 
igync) during which an electron loses half of its energy, would roughly be (in the observer's 
frame) tsync(-Ekev) = 1-2 x 10^ B"^/'^ 5~^/'^ ■, where E'kev is the observed energy in keV 
(Takahashi et al. 1996). The lag is defined as the difference of tsync(-^kev) from tsync(-E'o); 
which is approximately isync(-E^kev) when Eq is at much higher energies. Here we take Eq 
to be the logarithmic mean energy in the 3.0 — 7.5 keV band. Our subsequent analysis 
imphes 5 of 20 — 30, so we use the value of 25 as a plausible value (see § 3.2). We find the 
magnetic field B is 0.11 ± 0.01 G for 5 = 25. The solid line shows the model with the best 
fit parameter given above. 

2.2.2. Discrete Correlation Function 

We computed cross correlations using the discrete correlation function (hereafter DCF) 
of Edelson & Krolik (1988). We divided the 0.5 — 7.5 keV range into 5 energy bands and 
measured the time lag for each light curve compared to the 3.0 — 7.5 keV light curve. The 
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enor on the lag was determined from the uncertainty {la error) of the peak parameter 
obtained by the minimum fitting of the distribution of the DCF to a Gaussian. The 
results for both SIS and GIS data are shown in Figure 5 (d). Again the formula for the 
synchrotron coohng time was used to estimate the magnetic field strength (see § 2.2.1). 
The best fit value of B is 0.13 ± 0.01 G for 5 = 25 (sohd line). 

3. Time Dependent SSC Model 

3.1. Summary of the Numerical Code 

As we have seen in previous chapters, PKS 2155—304 showed rapid, large- amplitude 
variability throughout the ASCA energy bands. To describe the data in detail, it is essential 
to model the full time evolution of the synchrotron/Compton spectrum, incorporating 
radiative cooling processes, particle injection and escape. In the following, we consider a 
time-dependent one-zone SSC model, with a homogeneous emission region. Steady-state 
one-zone models have been shown to provide a good representation of the data for similar 
blazars Mrk 421 and Mrk 501 (Mastichiadis & Kirk 1997; Pian et al. 1998; Kataoka et al. 
1999a). 

Our kinetic code was developed as an application of the SSC model given in Inoue 
& Takahara (1996) and Kataoka et al. (1999a), which is accurate in both the Thomson 
and Klein-Nishina regimes. A spherical geometry with radius R is adopted. Our code is 
qualitatively similar to the one given in Chiaberge & Ghisellini (1999); in particular, we 
consider both the time evolution of the electron and photon distributions, as well as the 
effects introduced by the different photon-travel times in different parts of the emission 
region. For the steady state solution of the SSC code, we compared our results with Band & 
Grindlay (1985; 1986). The time evolution of the synchrotron spectrum was also compared 
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with the analytic solutions given by Kardashev (1962), for the simplest case of synchrotron 
cooling. The kinetic code will be described in detail in a forthcoming paper (Kataoka et al. 
1999b); here we briefly summarize the calculations, and present the results as applied to 
PKS 2155-304. 

We first consider a given electron population characterized by A^e(7,t), where 7 is the 
electron Lorentz factor and is the electron number density per unit volume per unit 
energy. The time t can be regarded as an initial time for the evolution of the electron and 
photon spectra. 

Given the electron population, the synchrotron emission including self-absorption can 

be calculated using the spherical solution to the radiative transfer equation, 

L,y..M = 47r^it:^^(l - —^[l - e-^ir. + 1)]) , (1) 

(Xi/ Til 

where jv and a^, are respectively the emission and absorption coefficients for synchrotron 
radiation (e.g., Blumenthal & Gould 1970; Gould 1979; Inoue & Takahara 1996). The 
optical depth in the plasma cloud along the central line of sight is = 2a,yR. 

We calculate the inverse Compton emission incorporating the effects of cross section 
reduction in the Klein-Nishina regime. The differential photon production rate q{e, t) is 

q{e,t) = j deon{eo,t) j d-fNe{-f,t)C{e,-f,eo) , (2) 

where eg and e are respectively the soft photon energy and the scattered photon energy 
in units of nieC^- The number density of soft photons per energy interval is n{eQ,t), and 
C(e, 7, eo) is the Compton kernel of Jones (1968). As the optical depth for Compton 
scattering is small, the self-Compton luminosity is obtained from the relation 

LsscM^fn'R'j!^^ , (3) 

where j^^^ is the emission coefficient of the Compton emission, related to q{e, t) as j^^*" = 
heq{e,t){A tt)"^ and h is the Planck constant. 



-13- 



The time evolution of the high-energy electrons in the magnetic field and the photon 
field is described by the following kinetic equation: 

qI ' = 5:^[(%ync + 7ssc)^^e(7, t)] + Q{l^ t) - , (4) 

where Q and ^esc are the injection rate and the escape time of the electrons, respectively. 
For simplicity, we set tesc to a constant value. The synchrotron and SSC cooling rates for a 
single electron are expressed as 

AaT'l^Us f f 

^'^"'^ " 3m c ' ^^^^ ^ J ^^^J ^^o^^^O' ^' ' 

where Ub is the magnetic field energy density. Thus the right-hand side of equation (4) is 
now described by the quantities at time t, and solved numerically to obtain the electron 
population at time i-|-At. We adopted an implicit difference scheme by Chang & Cooper 
(1970) to obtain non-negative and particle number conserving solutions. The iteration of 
the above prescription gives the electron and photon populations at an arbitrary time. 

Finally, we consider light-travel time effects. We proceed under the assumption 
that electrons are injected uniformly throughout a homogeneous emission region. Such a 
description may be adequate as long as the injection timescale is longer than R/c (Dermer 
1998), but is unphysical for shorter timescales since the particle injection process itself 
should take (at least) ^ R/c to inffuence the whole region. In a more realistic picture, a thin 
shock front may propagate through the emission region with a finite velocity Vs, supplying 
freshly accelerated electrons only in the front's vicinity (Kirk, Rieger & Mastichiadis 
1998). However, such a detailed calculation necessarily involves some additional, uncertain 
parameters, not to mention the assumption of a particular geometry. Instead, we will choose 
the injection to be uniform over the emission volume. Such a choice, albeit unrealistic, 
will allow us to clarify the important role of light travel time effects on blazar variability 
(We will, however, choose the duration of the injection to be comparable to R/c). The 
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observer would see, at any given time, photons produced in different parts of the source, 
characterized by an electron distribution of a different age. The observed spectra must then 
be a sum of the corresponding different spectra. Below, we will argue that light travel time 
effects may be essential in interpreting the shapes of blazar light curves. 

In order to incorporate the light-travel time effects properly in the calculation, we 
divide the source into {2tcis/ ^t) shces of AR thickness for each, where tcrs {— R/c) is the 
source light-crossing time and At is the time-step of the calculation. The interval At must 
be shorter than the shortest relevant time scale, e.g, the synchrotron or Compton cooling 
times or the injection time scale. We define AR = cAt. The schematic view of this division 
of the emission blob into shces is given in Figure 6. We first consider the slices in the source 
frame K', with line of sight perpendicular to the surface of the slices. In the observer's 
frame K, emission will be concentrated in the forward direction within a narrow cone of 
half-angle l/F (e.g., Rybicki and Lightman 1979), if the blob moves with a Lorentz factor 
r (~ 5 for our case). 

In our current assumption, the synchrotron cooling occurs homogeneously over the 
emission region, because the electrons are injected uniformly into the entire volume with a 
constant magnetic field strength. Even so, the Compton scattering should not be uniform, 
as each electron should experience the changes in the photon field after a different time 
interval depending on its position. The synchrotron photons should require a time ~ i?/c 
to be completely scattered, causing an additional delay for the response of the Compton 
photon spectrum. However, for this particular TeV blazar, synchrotron cooling must 
dominate over Compton coohng for the energy loss process of the electrons, as (1) the 
synchrotron luminosity is greater than the Compton luminosity, and (2) the reduction of 
the cross section in the Klein-Nishina regime significantly decreases the Compton scattering 
efficiency (see also § 3.3). In addition, the present gamma-ray data is not sampled well 
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enough to allow a comparison to the Compton model spectrum on very short timescales. 
Therefore, in calculating the Compton spectrum and energy loss, wc approximate that the 
synchrotron radiation (i.e., the soft photons for inverse Compton scattering) instantaneously 
fills the whole emission region, reducing the computational time (see also Chiaberge & 
Ghisellini 1999). 



3.2. Self-Consistent Solution for PKS 2155-304 

To reproduce the flare of PKS 2155-304 in 1994 May, we flrst model the multi- 
frequency spectrum from the radio to TeV bands, and determine the physical parameters 
relevant for the quiescent emission. Since the source activity is well-sampled, and appears 
to be rather smooth, we consider the level at ~ MJD 49492 as a locally quiescent state 
(sec Figure 1). We use an accurate solution of the one-zone homogeneous SSC model (§ 
3.1), to describe the multi-frequency spectrum of PKS 2155—304. To specify the spectral 
energy distribution for this source, we need the magnetic field B, region size R, beaming 
factor 5, escape time tesc, and the electron injection spectrum as input parameters. We 
adopted a specific form for the latter, Q{j) = Qe 7"^* exp(— 7/7inax), where 7inax is the 
maximum Lorentz factor of the electrons. For the minimum Lorentz factor of the electrons, 
7min, we set 7min = 1, although this choice is not important for our results. Thus seven 
free-parameters are required to specify the model. 

All these parameters are associated with seven observables: synchrotron peak frequency, 
Ps', Compton maximum frequency, Uc] synchrotron break frequency, i^br! variability time 
scale, ^var; synchrotron energy fiux, /sync; Compton energy fiux, /ssc! and radio (millimeter) 
spectral index, a. Despite the equal numbers of observables and parameters (7 observables 
for 7 parameters), the model cannot be specified uniquely, because the region size is 
described by an inequality R < ct^ax^. To determine model parameters uniquely, we make 
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the approximation it! ~ ctyg^T-S in the following. 

The parameters for the one-zone homogeneous SSC model are characterized by these 
observables as follows (see also Mastichiadis & Kirk 1997; Tavecchio, Maraschi & GhiseUini 
1998). The peak frequency of the synchrotron component in the observer's frame is given as 

1.2 X 10«S<57Lx • (6) 

We define the Compton maximum frequency Pc as the frequency where the Compton 
luminosity decreases by an order of magnitude from its peak observed flux (assumed here 
the highest energy bin in the EGKET data). For our PKS 2155—304 data, this corresponds 
to ~ 1 TeV and coincides with the VHE range (Chadwick et al. 1999). At this highest 
photon energy (~ hi/c), the emission from the TeV blazars are probably suppressed by the 
decreased Compton scattering cross section in the Klein-Nishina regime. Thus we obtain 
the approximate equality of the maximum electron energy and the photon energy 

Ima^rneC^S ~ hu^ . (7) 

The third relation is obtained from the ratio of the synchrotron luminosity to the 
Compton luminosity (see Kataoka et al. 1999a): 

^ 2 f 2 

R^cb'^ /BSC 

where dL is the luminosity distance. Combining equations (6) — (8), we can express i? as a 
function of only one parameter, 5. By equating R with ct^8, we can determine the relation 

between observables and 5, i?, B and 7max uniquely. 

To determine the escape time ^esc, another relation is required. An injection of a 
power law energy distribution of electrons up to a certain maximum energy (7max) into the 
radiating region should yield a steady-state electron distribution with a break in its index 
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at a characteristic energy (7br), determined by the balance between radiative coohng and 
advective escape (e.g., Inoue & Takahara. 1996). Thus at 7 = 7br, we expect tesc = -^cooi- 
This can be written as 

SnieC . . 

^^sc — 7777 ; 77 ^ ) (9) 

where Uson is the soft photon density, which is related to Ub via Us/Usoh = /sync//ssc 
in the Thomson regime. Thus we can express tesc as a function of B, S, and i/br- JA)r is 
determined as the frequency where the spectral index starts to deviate from that in the 
radio (millimeter) band, which presumably corresponds to the uncooled portion of the 
electron distribution. 

We summarize the 'input' observables in Table 1, and 'output' physical quantities in 
Table 2. The relation between the observables and model parameters is also shown in Table 
2. The remaining parameter Qe for injected electrons was determined to agree with the 
synchrotron luminosity /sync- The slope of the injected electrons s was determined simply 
by 2q; + 1. As a result, we obtain a magnetic field 5 = 0.14 G, region size R = 7.7 x 10^^ 
pc, and beaming factor S — 28. This result is quite similar to that reported in Tavecchio, 
Maraschi & Ghisellini (1998). They investigated the allowed parameter space for (B, 6), 
and conclude that B ~ 0.2 G with (5 ~ 25 is the most probable parameter set for PKS 
2155-304. 

To estimate the uncertainties of the parameters listed in Table 2 and the resultant 
model prediction of the multi-frequency spectrum, we considered different cases where the 
selected value of one observable is varied while the other six observables remain unchanged 
(see Table 1). We first changed the relation R — civar*^ to i? (1/3) ctva,i5. In this case, 
we obtain i? = 0.19 G and 6 = 37 with region size i? = 1.0 x 10^^ cm. Next we increased 
the Compton energy flux fssc by a factor of 3. This results in B — 0.12 G, S — 25 and R 
= 2.1 X 10^^ cm. Finally we increased the maximum Compton frequency by a factor of 2. 
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This changes the result most significantly, where B — 0.05 G, 5 — AO and R — 3.3 x 10^^ 
cm. However, when 5 ;^ 30, the TeV flux becomes considerably higher than that observed, 
because most of the Compton scattering takes place in the Thomson regime. Also note that 
any change of i4,r keeps all the parameters unchanged except for tesc- 

To obtain the steady state solution for electron and photon spectra based on the 
parameters hsted in Table 2, we start from the initial condition A^e(7, 0) = Qi'j) and 
calculate the time evolution of the spectra to more than t > 30 tcrs, assuming constant 
injection and escape. In Figure 7, we show the multi-frequency spectra of PKS 2155—304, 
as well as the calculated SSC model lines. The solid lines in Figure 7 show the synchrotron 
self-Compton spectra at the 'steady state', based on the parameters in Table 2. Although 
most of the data in Figure 7 were obtained non-simultaneously, variability is not large 
except for the X-ray band, so these spectra provide a reasonable constraint on the source 
parameters. One can see the discrepancy in the radio band, but at the higher frequencies, 
the model represents the observed spectrum quite well. The discrepancy in the radio 
band can be an effect of the emission arising from larger distances than the location of 
the keV/TeV emitting region. The one-zone models cannot account for this low-frequency 
emission, which is thought to be produced in a much larger region of the source (e.g., 
Marscher 1980). 

3.3. Modeling of the PKS 2155-304 Flare in 1994 

Based on the physical parameters selected above, and the resultant 'steady state' 
spectrum as an initial condition, we model the X-ray flare of PKS 2155—304 in 1994 May. 
Various types of flaring behavior were investigated by changing the parameters for the 
injected electron spectrum and/or the physical quantities in the emission region. In the 
following, we simulated the light curves, as well as the time evolution of spectra, when 
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7max increased by a factor of 1.6 during one t^rs interval from the start of the flare. After 
the calculation was performed in the source frame, it was transformed into the observer's 
frame for comparison with the observational data. We take the time-step to be At = 2000 
sec throughout the calculation. This corresponds to At < 100 sec in the observer's frame, 
much shorter than the variation time scales, such as the synchrotron cooling time (see § 2). 

We varied the injected electron spectrum as Qei'j) — Qe exp(— 7/1.67inax) for < 
t < tcrs and Qeil) — Qe 6xp(— 7/7inax) Otherwise. In Figure 8, we show the calculated 
light curves from EUVE to ASCA energies. The flux was normalized to that for the steady 
state {t — 0), and the time axis was normalized to the source light-crossing timescale {R/c). 
The symmetric light curve during the flare is reproduced quite well. Notably, the peaking 
time of the flare at lower energies lags behind that for higher energies, and the amplitude 
of the flare becomes larger as the photon energy increases, which agree qualitatively with 
the observational data. 

The position of the peak time is determined by the balance of slices in which the 
emitted flux is increasing and slices in which the decaying phase has already started. For 
example, at the highest X-ray energy band which corresponds to 7 ~ 7max, the peak of the 
light curve will occur at a time t, where tcrs < t < 2tcrs- This is because these electrons in 
the slices cool much faster than tcrs and the volume which contains the flare information 
is maximum at the center of the sphere, which can be observed after tcrs from the start of 
the flare (see Figure 6). The whole emission region requires 2 tcrs to be completely visible 
to the observer. This will cause a decaying of the flux (foreside slices) and an increasing of 
the flux (backside slices) at the same time. However, in the lower energy band where tcooi 
^ tcrs, the electrons do not cool effectively. The emitted flux will continue to increase even 
after 2 tcrs- This combination of the increasing/decreasing phase of slices can cause a time 
lag in the position of the peak as can be seen in Figure 8. 
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The dashed lines in Figure 2 (a) and (b) show the spectral evolution calculated from 
the adopted model. The 2 — 10 keV flux was obtained by integration of the calculated 
spectrum in the 2 — 10 keV band, so it can be readily compared with the observational 
data. The photon spectral index was simply determined from the ratio of the fluxes at 0.7 
keV to 7.5 keV. Since the ASCA observation started after the onset of the flare, we shifted 
our simulated light curves in the time-axis for comparison. The fact that the peaking time 
of the photon index leads that for the flux is quite well reproduced. This can be also seen 
in Figure 3, where the 'clockwise' hysteresis in the flux (2 — 10 keV) versus photon index 
(0.7 — 7.5 keV) plane can be clearly seen. 

To make a quantitative comparison of the observed and modeled light curves in 
different energy bands, we analyzed the light curves in the same way as that for the 
observational data. The dashed lines in Figure 5 (a) — (c) were calculated from a Gaussian 
flt of the simulated light curves, converted to the observer's frame. Figure 5 (a) shows the 
duration of the flare, determined from the standard deviation (a) of the Gaussian. This 
stays nearly constant at ~ 1.5 x 10^ sec, but the model shows a sign of increase at lower 
X-ray energies - a for 0.5 keV is 8 % longer than that for 5.0 keV. Although we do not 
detect any such increase (broadening) in the observed X-ray light curves (Figure 1), flares 
observed by EUVE and lUE showed considerably longer time scales than that in the ASCA 
band (Urry et al. 1997). This may suggest that the duration of the flare actually increases 
at lower energies, if both flares have the same origin. Also note that the Gaussian flt of 
the simulated light curves could make systematic errors on a, if the symmetry of the light 
curves is broken only in the lowest X-ray energy bands (see also the discussion in § 4). 

The amplitude of the flare {Ap = Ci/Cq) was calculated to be 0.5 for 0.5 keV and 1.6 

for 5 keV - precisely in agreement with the observational data (Figure 5 (b)). Wc compared 
the light curves at 5 keV to those for the lower energies, and found a 'soft-lag' where the 0.5 
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— 1 keV photons lag behind 5 keV photons about 4 ksec, quantitatively in agreement with 
the observational data (Figure 5 (c)). Importantly, this result implies that the observed 
time lag is well represented by the difference of the synchrotron cooling time scales. 

We also computed the discrete correlation function for the model light curves. The 
result is given in Figure 5 (d). Again, we obtained a similar result as with Gaussian fits, 
and verified a 'soft-lag' where the 0.5 — 1 keV photons lag behind 5 keV photons about ~ 
4 ksec. 

Finally, in Figure 7 we show the time evolution of the multiwavelength synchrotron 
self-Compton spectra after the start of the flare, at t/tcrs — 1-0, 1.6 and 2.4 (for comparison, 
see also Figure 8). The X-ray spectrum clearly becomes harder when the flux increases. 
Also note that the flux variation of the Compton spectrum is smaller than that for 
synchrotron spectrum. This is because the reduction of Compton cross section in the 
Klein-Nishina regime strongly suppresses the increase of the flux at the highest photon 
energies (~ hUf.). 

We also investigated other scenarios for the flare, such as increasing the normalization 
of the injected electrons (^e) or the magnetic fleld B strength, but these did not give good 
representations of the data. In a more realistic situation, all of these parameters probably 
vary simultaneously, but our modeling implies that the most essential parameter during the 
flare is 7max- We will discuss more about this in the next chapter. 

4. Discussion 

Our calculation shows that a one-zone homogeneous SSC model incorporating time 
evolution of the electron and photon distributions can be a viable tool to study the rapid 
variability of blazars. Investigating both the light curves and the spectral evolution of the 
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X-ray flare of PKS 2155—304 in 1994, our modeling reproduces the observational data quite 
well, despite the relatively simple assumptions regarding the flare mechanism. 

All the physical parameters for the model were self-consistently determined from seven 
observables (Table 1). We found that the magnetic fleld B = 0.14 G, region size R = 7.7 
X 10~^ pc, and beaming factor S = 28 give a good representation of the observational data 
(Table 2 and Figure 7). Although the present spectrum is too sparsely sampled, especially 
in the hard X-ray and 7-ray bands, we showed that our derived parameters are not affected 
by more than a factor of 2 even if taking the uncertainty of the observables into account. 

The magnetic field derived from the multiband spectrum is nicely consistent with that 
calculated from the X-ray variability based on the synchrotron cooling model (§ 2); while 
the Gaussian (§ 2.2.1) and DCF (§ 2.2.2) fits give slightly different values for B, both 
yielded the same value within the uncertainties, S ~ 0.1 G. Regarding this point, we should 
note that the errors on the lag determined from the DCF are somewhat ambiguous, because 
they are calculated from the errors on the peak of the distribution of DCF. When the DCF 
distribution is not well represented by a Gaussian, this can cause additional errors in B, 
which may account for the difference in the magnetic field determined from the Gaussian 
fit of the light curves. 

We also found that in the X-ray bands, the duration of the fiare is expected to be 
nearly constant in energy, forming very symmetric time-profiles. We considered a model 
where rapid local variability is relaxed (smoothed) on a light crossing time scale (e.g., 
Chiaberge & GhiseUini 1999; Giommi et al. 1998). In this scenario, the rising and falling 
phases of the flare are characterized by the same time scale R/c, forming a quasi-symmetric 
shape of the light curves. This can be clearly seen in Figure 1 and Figure 4, where the 
observed fight curves are well represented by a simple Gaussian function (see § 2.2.1). The 
standard deviation a determined from the Gaussian fit (see Figure 5 (a)) is thus expected 
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to be on the order of ~ R/c5, in the observer's frame. Importantly, the synchrotron coohng 
time is calculated to be ~ 7 x 10^ sec for 1 keV photons (see § 2.1.1) which is shorter than 
a (assuming B = 0.11 G and 6 = 25). However, at lower energies where tgync > R/cS, it is 
possible that the symmetry of the light curves is broken and a longer time scale decay could 
be observed - this is why the model hne in Figure 5 (a) begins to rise in the lowest X-ray 
energy bands (^ 1 keV). 

Effects other than that due to hght travel time could also make the hght curves 
symmetric, having constant duration independent of the photon energy. For example, 
Dermer (1998) considers a similar but simpler model, assuming that the electron injection 
event continues for at least a light-crossing time. In describing the time evolution of the 
electron/photon spectra, he discusses quantities which are averaged over timescales shorter 
than R/c. Highlighting the case of an s = 2 electron injection power- law index, he found 
that the light curves can become symmetric, with energy-dependent lag in the mean time 
of the flare in various energy bands. Importantly, he showed the flare mean time stays 
nearly constant in energy if tcooi is less than or comparable to the duration of injection, as 
is the case for our modeling. However, it is to be noted that symmetric light curves can be 
reproduced in his model only when the particular electron injection distribution of s = 2 is 
selected (see equation (2) of Dermer 1998). 

We believe that relaxing of rapid, local variabihty by light-travel time effects must 
generally be important in blazar variability. As noted above, if the only timescales relevant 
in the light curves are those associated with cooling and injection, we should mainly 
observe asymmetric light curves, except for the special case of s = 2. However, flares with 
symmetric time-proflles were observed many times from this particular source previously, 
in various states of activity with X-ray differential photon indices ranging from 2.0 to 
3.5 (e.g., Sembay et al. 1993; Edelson el al.l995; Giommi et al. 1998; Chiappetti et al. 
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1999), perhaps reflecting a range of electron injection indices. To account for the general 
variabihty trends of this source, it thus seems preferable to include light-travel time effects. 
Also note that symmetric light curves are not unique phenomena for PKS 2155—304, but 
also commonly found for other TeV blazars (e.g., Takahashi et al. (1999) for Mrk 421). 

To describe the PKS 2155—304 flare in 1994 May, we adopted a model where the 
flare is due to an increase of 7inax which continues for one source light-crossing time scale, 
^inj = ^crs- Regarding this point, we note that while longer durations of electron injection, 
i.e., tinj ^ tcrs, may be physically possible, this cannot have been the case for our PKS 
2155—304 flare. When the injection is much longer than tcrs, both the electron and photon 
distributions will have enough time to reach a 'new' equilibrium state, forming a plateau 
in the light curve (e.g., Mastichiadis & Kirk 1997; Chiaberge & Ghisellini 1999). Such a 
plateau, although not as common, has occasionally been observed in the case of Mrk 421 
(see, Takahashi et al. 1999; Kataoka et al. 1999b for more detailed discussion), but not in 
our PKS 2155—304 data. If the injection results from a shock propagating with velocity 
Vs through the emission region, the duration of the injection event should be characterized 
by the shock crossing timescale R/vs (e.g.. Kirk, Rieger & Mastichiadis (1998)). Thus our 
current assumption of t^j — tcrs corresponds to a highly relativistic shock, Vs ~ c. 

Despite the simple description of the injection (an increase of 7inax with other 
parameters fixed), it provides a good representation of the observational data. More 
realistically, this may correspond to a new population of electrons injected with a harder 
distribution during the flare. In the picture of shock acceleration, electrons with larger 7 
require a longer time to be accelerated (e.g.. Kirk, Rieger, & Mastichiadis. 1998), so our 
assumption of an increase in 7niax can be a good approximation as long as the acceleration 
timescale tacc at jmax is well shorter than tcooi- 

Our numerical calculation also predicts the flare amplitude, duration, and the time 



-25- 



lag of the same flaring event in the lUE and EUVE bands. The agreement with the data 
(Urry et al. 1997) is quahtatively correct but quantitatively there are discrepancies of a 
factor of ~4 if we interpret these flares as having the same origin. The time lags from the 
highest ASCA band (~ 7 keV) are estimated to be AtEuvE — 0.2 day and A.tiuE — 0.5 day, 
compared to the observed values of ~1 day and ~ 2 days, respectively. However, the EUVE 
'peak' is very flat, and both EUVE and lUE flares could conceivably represent multiple 
injections. The flare amplitude is estimated from the model to be Ap (= Ci/Cq; see § 2.2.1) 
= 0.2 for EUVE (compared to ~ 0.5 observed) and Ap = 0.05 for lUE (compared to ~ 0.35 
observed). The duration of the flare is calculated to be 0.2 day for EUVE and 0.5 day 
for lUE . These values are nearly consistent with EUVE observational data, but disagree 
more strongly with lUE data. These discrepancies may indicate that another mechanism 
is required to account for the low-energy flare, or that there are multiple smaller flares 
uncorrelated with X-ray variations, or that the one-zone homogeneous assumption can not 
be applied at lower energies. 

Recent work by Georganopoulos & Marscher (1999) also flnd it difficult to reproduce 
the flare amplitude observed with lUE in 1994 May, even assuming a more realistic, 
inhomogeneous SSC jet model with many parameters. They suggest that one of the origins 
of the discrepancy could be a mild re- acceleration of the electrons in the injected plasma. 
Such a mild re-acceleration may be manifested mainly at lower energies, thereby increasing 
the flare amplitude relative to that at higher energies. However, the short coverage of the 
lUE , EUVE and ASCA data during the 1994 May campaign is insufficient to deduce 
further conclusions. 
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5. Conclusions 

The TeV blazar PKS 2155-304 was observed with the X-ray sateUite ASCA in 1994 
May, as part of a multiwavelength campaign. A rapid, large- amphtude flare was detected at 
the beginning of the observation. During the flare, the change in the hard X-ray flux led the 
change in the soft X-ray flux by ~ 4 ksec, as derived by two different methods. We showed 
that the hght curves can be fltted with a Gaussian plus constant offset. The magnetic held 
B was estimated to be ~ 0.1 G for 6 = 25, assuming the energy-dependent delays are due to 
the different synchrotron lifetimes of the relativistic electrons. We also found that the flare 
amplitude increased as the photon energy increased, while the duration was nearly constant 
throughout the ASCA band. Using a time dependent model of synchrotron self-Compton 
emission in a homogeneous region, we considered a scenario where the electrons are 
uniformly injected into the emission region, and rapid, local variability is smoothed out 
by light-travel time effects. In our modeling, the lag of the peaking time can be explained 
by the different combination of the increasing/ decreasing phase of local emission regions 
("slices"). The spectral energy distribution of PKS 2155—304 from the radio to TeV bands 
was well represented by the parameters S ~ 0.1 — 0.2 G, i? ~ 10"^ pc and 5 ~ 20 — 30. 
Importantly, these parameters were self-consistently determined from seven observables. 
We found that the X-ray flare in 1994 May is well explained by a change in the maximum 
energy of the electron Lorentz factor (7max) by a factor of 1.6. However, the relation to the 
lUE and EUVE flares remains uncertain because of the short coverage of the observations. 
More data, especially longer and continuous coverage at neighboring frequencies such as 
UV to X-rays, are necessary to fully understand this source, and the physical processes 
operating in blazars. 
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Fig. 1. — Time history of the X-ray emission of PKS 2155-304 during the 1994 May 
campaign. Upper Panel: Combined SIS hght curves in the lower X-ray band (0.5 — 1.5 
keV). Middle Panel : Combined SIS light curves in the the higher X-ray band (1.5 — 7.5 
keV) . Both light curves are normalized to their average count rate. The arrows indicate the 
time interval used for the Gaussian fit and discrete correlation function (DCF) in § 2. The 
best fit Gaussian is superposed on each curve to guide the eye. Lower Panel: Time history 
of the hardness ratio, defined as the ratio of the SIS count rates at 1.5 — 7.5 keV to those 
at 0.7 - 1.5 keV. 

Fig. 2. — Detailed time history of PKS 2155-304 during the 1994 May campaign. Each 
data point corresponds to an equal 5 ksec interval and all SIS/GIS data are combined for 
the fit. The model is a power law with free absorption, (a): Variation of the differential 
photon index in the 0.7 — 7.5 keV band. The dashed line is a model prediction as described 
in § 3. (b): Variation of the 2-10 keV flux in units of 10~^° erg cm"^ s~^. The dashed line 
is a model prediction discussed from § 3. 

Fig. 3. — Evolution of the X-ray spectrum of PKS 2155—304, in the flux versus photon 
index plane. Arrows indicate the evolution during 1994 May observation. A 'clockwise loop' 
is clearly seen. The solid line connects the observational data, while the dashed line is the 
model prediction given in § 3. 

Fig. 4.— A Gaussian flt to the time proflle of PKS 2155-304 during the flare (May 19.2 
— 20.0 UT; see also Figure 1). The count rates of both SIS detectors are summed. The 
solid line corresponds to the best fit with f{t) = Cq + Ci x cxp {-{t - tp)V2a2), while the 
dashed line is the constant offset Co (see § 2.2.2). Upper Panel: Light curve in the 0.5 — 1.0 
keV band; Xred — ^-^ ^9 dof. Lower Panel: Light curve in the 2.0 — 3.0 keV band; x^^^j 
= 1.1 for 49 dof. 



-32- 



Fig. 5. — Measurement of the source parameters for a homogeneous synchrotron model 
describing the 1994 May flare of PKS 2155—304. The circles are combined SIS data, while 
the crosses are combined GIS data, (a): The flare duration, estimated by the standard 
deviation, cr, of the Gaussian fits at various X-ray energies (see § 2.2.1). The sohd line is a 
constant fit (cr = 1.5 x 10^ sec), while the dashed line is calculated from the Gaussian fit 
of the model light curves in § 3. (b): The flare amplitude {Ap = Ci/Cq) at various X-ray 
energies, determined by a Gaussian flt in § 2.2.1. The dashed line is calculated from the 
Gaussian fit of the modeled hght curves in § 3. (c): Time lag of photons of various X-ray 
energies versus the 3.0 — 7.5 keV band photons, calculated from the peaking time determined 
by a Gaussian fit. The solid line is the best fit with a function given in § 2.2.1 (5 = 0.11 G, 
5 — 25). The dashed line is calculated from the Gaussian fit of the model hght curves in § 3. 
(d): Time lag of photons of various X-ray energies versus the 3.0 — 7.5 keV band photons, 
calculated using the discrete correlation function (DCF). The solid line is the best fit with a 
function given in § 2.2.1 [B — 0.13 G, 5 — 25). The dashed line is as in (c). 

Fig. 6. — The schematic view of the division of the emission blob into 'slices'. Left Panel: 
Source frame {K'). Right Panel: Observer's frame {K). The emission blob is first cut into 
slices in the source frame and contributions from each shell are summed. In the observer's 
frame, the radiation is concentrated in a narrow cone with a half angle 9 ~ l/F. 

Fig. 7. — Multi-band spectrum of PKS 2155-304. The filled circles are the nearly 
simultaneous lUE and EUVE data reported in Urry et al. (1997). The filled squares are 
the X-ray data obtained with ASCA (this work). The data from both GISs are combined. 
The high state ASCA data are from the 15 ksec of the peak fiux of the fiare, while the low 
state data are from the 15 ksec near the end of the observation. EGRET data and TeV data 
are from Vestrand, Stacy, & Sreekumar (1995) and Chadwick et al. (1999), respectively. 
The other non-simultaneous data are from the NED data base. The sohd fine and dashed 
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lines represent the time evolution of the photon spectra calculated from the time-dependent 
SSC model described in § 3. The physical (model) parameters are given in Table 2. The 
multi-frequency spectra at t/tcrs = 0, 1.0, 1.6 and 2.4 are shown respectively. 

Fig. 8. — Upper Panel: Simulated light curves at different UV/X-ray energies, reproducing 
the rapid flare of PKS 2155—304 in 1994 May. We plotted the time evolution of the fluxes at 
different frequencies, normalized to the initial {t — 0) value. The quasi-symmetric shape of 
the high energy light curves and the increasing time lag of the peak with decreasing energy 
is clearly seen. Lower Panel: Assumed change of 7niax during the flare. The 7max is assumed 
to increase by a factor of 1.6 from the initial value for < i < icrs- Time-axes of both panels 
are in units of the source light-crossing time (icrs)- In the observer's frame, tcrs corresponds 
to 0.3 day. 
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Table 1. Input Observables for PKS 2155-304 



Parameters 




Input Values 


z/s,i7 (in 10^^ Hz) 




0.3 


i^c,27 (in 10^'' Hz) 




0.3 


Z4,r,i4 (in 10^^ Hz) 




0.5 


is (in 10"'' sec) 




0.3 


/sync (in 10"^" erg cm"^ 


s-i) 


13 


/ssc (in 10"^° erg cm'^ 


s-i) 


4.1 


a 




0.18 



Input parameters (observational quantities) 
determined from the multi-frequency spectrum 
of PKS 2155-304 (Figure 7), as defined in §3.2. 
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Table 2. Output Physical Quantities for PKS 2155-304 



Parameters 


Relation to Observables 


Output Values 


S 

B{G) 
R (cm) 

^esc (sec) 
'Ymax 


3 0x10-6//* fl/i ---1/8 -1/4 1/2 -1/4 
O.UXiU /sync/ssc '^8,17 '^c,27'-5 

3 8x10-9//"^^/" .-1/8 3/4 -3/2 -1/4 
O.OXIU /sync/ssc '^s,17'^c,27 ''5 

S Svin9//4fV4 ,-1/8 -1/4 1/2 3/4 
S.SXiU /sync/ssc ^s,17 ^c,27^5 

6.3xlO"rf-VVs;ycVs'^c(l+/sSc//sync)-^<f/VS</f4/^ 

2 7X 1012^-^4 ,1/8 1/4 1/2 1/4 
Z./XIU Uj^ /sync /sSC'^s,17'^c,27''5 


28 
0.14 
2.4x10^6 
8.6x10^ 
8.3x10^ 


Qe (cm~^ S~^) 


— (normalized to agree with /sync) 


6.1x10"^ 


S 


2a+l 


1.35 



Output model parameters for the multi- frequency spectrum of PKS 2155—304 (Figure 7). (Il 
= 1.47x10^'' cm is the luminosity distance of PKS 2155—304. The emission region size R is 
approximated as ~ ctyaiS. We adopted a cutoff power law for the injected electron population : 
Qil) = Qe exp(— 7/7max)- We set the minimum Lorentz factor of the electrons, 7inin = 1, for 
simplicity. 
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